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1 Introduction 

Let M be a Poisson manifold endowed with a bivector Po fulfilling the Jacobi condition 

[Po,Po] =0 

with respect to the Schouten bracket [. , .] [15] . 

A bi-Hamiltonian manifold M is a smooth (or complex) manifold endowed with two com- 
patible bivectors Po , Pi such that 

[Po,Po] = [Po,Pi] = [Pi,Pi] = 0. (1.1) 

Classification of compatible Poisson bivectors and the corresponding bi-integrable systems with 
integrals of motion 

{Hi,Hj} = {Hi^^x = 0, i,j = 1, ...,n, 

is nowadays a subject of intense research [3] [13l [14] . Of course, direct solution of the equations 
([Lip is generally quite difficult. We can try to lighten this work by using properties of the given 
Poisson manifold (M, Po). 

Bivectors Pi fulfilling the compatibility condition [Po,Pi] = are called 2-cocycles in the 
Poisson-Lichnerowicz cohomology defined by Po on M [6]. The Lie derivative of Po along any 
vector field X on M 

P 1 =C X (P ) (1.2) 

is 2-coboundary, more precisely it is 2-cocycle associated with the Liouville vector field X. For 
such bivectors Pi the system of equations (jl.ip is reduced to the single equation 

[£ x (P ),Cx(Po)] = 0, & [£|(P ),P ] = 0. (1.3) 

The second Poisson-Lichnerowicz cohomology group Hp (M) of M is precisely the set of 
bivectors Pi solving [P ,Pi] = modulo the solutions of the form Pi = Cx{Pq)- We can 
interpret Hp (M) as the space of infinitesimal deformations of the Poisson structure modulo 
trivial deformations. For regular Poisson manifolds this cohomology reflects the topology of the 
leaf space and the variation in the symplectic structure as one passes from one leaf to another 

The aim of this paper is to show some particular solutions of the equation (|1.3p and to 
discuss separation of variables for the corresponding bi-integrable system on so* (4) . 
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2 The Lie-Poisson bivectors on so* (4). 

Let us consider the semisimple Lie algebra so(4). The dual space M = so* (4) is the Poisson 
manifold [15 . Since M — so* (4) = so* (3) ©so* (3), we can introduce the following coordinates 
z = (s.t) on M, where s = (sx,S2,S3) and t = (tx,t 2 ,t 3 ) are two vectors in R 3 . As usual we 
identify (K 3 , A) and (so(3), [.,.]) by using the well known isomorphism of the Lie algebras 

/ o Z 3 -22 \ 

Z = [Zl, Z 2 , Z 3 ) — > Z M = -^3 zi , (2.1) 
\ z 2 —z\ / 

where A is the cross product in R, 3 and [., .] is the matrix commutator in so(3). In these 
coordinates the canonical Poisson bivector on so* (4) is equal to 

The generic symplectic leaves are the level sets of two globally defined Casimir functions 

3 

C x = (8,8) = |s| 2 c 2 = (t,t), P dC h2 = 0. (2.3) 

i=X 

In the paper [TO] the overdetermined system of equations (jl.lj) on so* (4) has been directly 
solved in the class of the linear Lie-Poisson bivectors Pi, i.e using the nzats P\° = ^ a m z m for 
the components of the second bivector. 

According to [TU] the Frahm-Schottky, Steklov and Poincare systems are bi-integrable 
systems and the corresponding second bivectors Pi — Cx(Po) are generated by the Liouville 
vector fields Xj — X^di with the following components 

Xi = (oi*i, a 2 t 2 , a 3 t 3 , aisi, a 2 s 2 , a 3 s 3 ), 

X 2 = (Xi,aX 2 ) a, a\, a 2 , a 3 G C, 

X 3 = (0,0,at 3 ,0,0,-as 3 ), 

where 

Xx = iaitl + ^H Sl ,a 2 t 2 + ^±^s 2 , a 3 t 3 + ^|±^s 3 ), 

X 2 = (-a 2 a 3 S! - ^^-tx, -axa 3 s 2 - ^^t 2 , -axa 2 s 3 - ^^t 3 ), 
The remaining four solutions of the system are associated with the following vector fields 

Xi = {{a 3 + a 2 )sx, (a 3 + ai)s 2 , (ax + a 2 )s 3 , (b 2 + 63)^1, {bx + h)t 2 , (bx + b 2 )t 3 ) 

X 5 = (-fsi + Mi, -f s 2 + b 2 t 2 ,b 3 t 3 ,bxsx - pi,b 2 s 2 - ft 2 , (b 3 - a)s 3 - at 3 ) 

X 6 = (a 2 sx,axs 2 , (ax + a 2 )s 3 - bt 3 , ctx, ct 2 ,0) 

X 7 = (asx,as 2 , 0, cis 3 + bti, c 2 s 3 + bt 2 , 0), 

where bx — ±b 2 for the vector field X§. The corresponding bi-integrable systems with quadratic 
integrals of motion have been completely described in 10J . 

We have to keep in the mind that the vector fields Xk and the corresponding bivectors 
Pi = Cx k Po are defined up to canonical transformations of the vectors s and t, which preserve 
the form of the canonical bivector Pq (12. 2p [5] . 
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3 Quadratic Poisson bivectors on so* (4) and e*(3) 

After the linear Poisson structures, it is natural to look at quadratic structures. Substituting 
the following anzats for the components of the Liouville vector fields X m = Y^>j a m z i z j m t° 
(|1.3|) one gets a highly overdetermined system of quadratic equations on the 126 coefficients 
a\ J n . Unfortunately, we cannot get and classify all the solutions of this system even by using 
modern computers and modern software. 
So, in this section we will suppose that 

P 1 dCi ;2 =0. (3.1) 

It means that symplectic leaf of P\ are contained in those of Pq. In general bivector Pi could 
have some more Casimirs, so that their symplectic leaf could be smaller. For example, if X = 0, 
then equation (|3.1[) is true, but the symplectic foliations of Po and Pi are different. However 
in this case Pi = does not linear or quadratic bivector on so* (4). 

Proposition 1 For all the quadratic bivectors Pi — Cx{Po) on so* (4) restriction i3.1\) leads 
to quadratic Poisson structures having the same foliation by symplectic leaves as Pq only. 

We will prove this Proposition by using so-called Darboux-Nijenhuis variables. 

The additional restriction on Pi = Cx(Po) is a linear equation with respect to X, which 
allows us to get solutions of the following system of equations 

Cx(Po)dCi,2 = -PodX(C 1 , 2 )=0 and [C X {P ), Cx(P )] = (3.2) 

Any solution X of this system (|3.2|) gives rise to a Poisson bivector Pi = Cx(Po) fulfilling 
equations and (|3.ip . 

We solved equations (|3.2p in the class of quadratic vector fields X = X m d m with the 
components X m = Yli>j a m z i z j by using one of the modern computer algebra system and got 
the following three solutions. 

Proposition 2 Let a and b are numerical vectors and a, (3 are arbitrary parameters. The 
vector field X — (Xi,X 2 ) with components 

Xi =a(sA(aAs)) + 0{b, t)(s A a), X 2 = a(t A (6 At)) 

generates a first bivector fulfilling equations and i3.1\) 

(i) = ( 2a(a, S )s M /3[(aAs)® (6At)]\ 

1 U[(aAs)®(fcAt)] T 2a{b,t)t M )' K '' 

Using orthogonal transformations of the vectors 

s -> s' = Uis, and t -> t' = U 2 t (3.4) 
where Ui j2 are orthogonal matrices, we can always put 

a = (0,0,03), b= (0,0, b 3 ). 

Using scaling transformation 

P^^XP^, AeM, (3.5) 

we can put a = 1 or /3 = 1. 

Proposition 3 Let a and b are two complex vectors, such that 

(a,b) = (b,b} = 0, 
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where (x,y) — Y^—i x iUi * s the inner product of two vectors. The vector field X — (Xi,X%) 
with components 

X\ = ~ s A (a A s), 

X 2 = -\t A (a At) +i\a\~ 1 (a,s)(a At) - i{b, s)(b At), 
generates a second bivector fulfilling equations and i3.1\) 

p (2) _ ( (a,s)s M -i[\a\- 1 (aAs)(»(aAt)-(bAs)®(bAt)]\ 

\il|a| 1 (aAs)®(aAt)-(bAs)®(bAt) I — (a,t)t M J 



here \a\ = \/ a\ + a\ + a| and (x (8> y)ij = a^yj- 

Using orthogonal transformations (13. 4p and scaling (|3.5|) we can always put 

a = (0,0,1), 6= (6i,«6i,0). 

Proposition 4 Let a, 6 and c are t/iree complex vectors, such that 

(a, 6) = (a, c) = (6, 6) = (c, c) = 0, b A c ^ 0. 
T/ie vector field X = (X\,X2) with components 
X\ — |sA(aAs), 

= ± t A (a A t) - i|a| _1 (a, s)(a At) - s)(c A i), 
generates a third bivector fulfilling equations il.l)) and iS.l]) 

(3) _ ( {a,s)s M i[\a\- 1 (aAs)®(aAt) + (bAs)®(cAt)] \ 

\—i\\a\ (aAs)(g>(aAt)+(6As)®(cAt)J {a,t)t M J 

Using orthogonal transformations (|3.4[) and scaling ()3 . 5|) we can always put 
a =(0,0,1), b= (bi,ibi,0), c = (ci, -ici,0) . 
Let us consider the canonical bivector on M = e*(3) 

p =L° r) (3.8) 



Xm Jm 

and the corresponding Casimir functions 

Ci = |x| 2 — x\+ x\ + x\, C 2 — (x, J) = x\ Ji + X2J2 + X3J3. 

On M = e*(3) the system of equations (|3.2p has only one nontrivial solution in the class 
of quadratic vector fields X = ^2 X rn d m with the components X m = 53i>j ; a m z i z .j? where 
z = (a;i,x 2 ,a;3, Ji, J 2 , J 3 )- 

Proposition 5 If a and b are two vectors, such that \a\ = 0, and a is arbitrary parameter, 
then the following bivector on M = e*(3) 

(a.x). (xAJ)®a+{a.J)x M + -2( — +oi{a.b))(x®x—\x\ 2 ')+a(bAx)®(aAx) 

P/ 43 = [ 1 (3.9) 

{a.J).J M + (b,x)x M + ^(^+a(a.b)){xA.J) M — ( (a Arc) A (b A J) J 



satisfies equations 11. 1]) and liS.l]) 

This bivector is a Lie derivat 
the corresponding Liouville vector field X. 

Poisso 

(2) (3) 

foliation with P and that the Poisson bivectors P x and P± are equivalent 



This bivector P^ is a Lie derivative of Pq . For the brevity we omit the explicit expression for 
In the next Section we prove that the Poisson bivectors P^ 171 ' have a common symplectic 
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4 The Darboux-Nijenhuis variables 



Let us consider a bi-hamiltonian manifold M with the non degenerate Poisson bivectors Pq and 
Pi. By definition a set of local coordinates (qi,Pi) on M is called a set of Darboux-Nijenhuis 
coordinates if they are canonical with respect to the symplectic form 

n 

lu = Pq 1 = dpi A dqi 

i=l 

and put the recursion operator N = PiP^ 1 in diagonal form, 



„ ^ / d J d J \ 



(4.1) 



This means that the only nonzero Poisson brackets are 

{?i,Pj}o = <%, = q&j, {Qi,Qj}o,i = {PhPj}o,i = 0. (4.2) 



According to (|4.1[) coordinates <?i are eigenvalues of N, i.e they are roots of the minimal poly- 
nomial .4(A) = (det(JV - AI)) 1/2 of N. 

In order to get momenta pi^ 2 we can directly solve equations (|4.2|) with respect to the 
functions pi$(s, t). The standard computation problem is that variables p\$ are defined up to 
canonical transformation pi —> Pi + fi(qi), where ft are arbitrary functions on qi. 

In order to construct recursion operator N on the generic symplectic leaves of so* (4) we 
will use analog of the Andoyer variables [TJ [5] • 

4.1 Analog of the Andoyer variables on so* (4) 

Let us introduce the following analog of the Andoyer variables [L 

,4 i I s 2 + isi+t 2 + iti \ 

ui = s 3 + £3, Vi = —i In — . (4.3) 

\^{ti + siy + { S2 + t 2 y) 

and 



u 2 = y/d + C 2 + 2s x ti + 2s 2 t 2 + 2s 3 t 3 , 

^3^*1 - S 3 C 2 + (*3 - S 3 )(siti + S 2 t 2 + S 3 t 3 ) 



(4.4) 



v 2 = arccos 



\/(h +si) 2 + (s 2 +t 2 ) 2 y/CiC 2 - (sih + s 2 t 2 + s 3 t 3 ) 2 J 
The inverse transformation in coordinates 

Ji = Si+U, Xi = x(si—ti), xeC, (4.5) 

looks like 
and 



Ji = \ u 2 — u\ sin vi, J 2 — \ u 2 — u( cos vi , J 3 = u 



Xi = — 2 1 - smui + \ L 2 — x z u 2 j- smv 2 cosui 



x 2 = — ! ! cos vi + \/C 2 — x 2 u 2 — % sin v 2 sin vi 



x 3 =M - a/1 - £*/C 2 - _ §cos V2 
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Here 



Ci = x\ Ji + x 2 J 2 + x 3 J 3 = k(Cx - C 2 ), 



C 2 = x 2 +x 2 +x 2 + x 2 (J 2 + J 2 + J 2 ) = 2>c 2 {C 1 +C 2 ). 
In (a;, J)-variables (|4. 5[) canonical bivector Pq (|2.2p on so* (4) reads as 

p _ f x 2 J M xm 

° V XM Jm 

At x — * this bivector is reduced to the canonical bivector Pq on e*(3), which is a dual space to 
the algebra e(3) of Euclidean group P(3). After contraction x — > our variables u, v (j4.3H4.4p 
coincide with the Andoyer variables on e*(3) PQ. 

The projection of the canonical bivector P on the generic symplectic leaves of M — so* (4) 
or M = e*(3) in (u, u)-variables (|4.3H4.4p looks like 

/ 1 \ 
1 
-10 
\ -10 0/ 

Using this analog of the Andoyer variables we can easy obtain the projection Pi of any bivector 
Pi fulfilling (I3~T1) . 

So, we can introduce symplectic form u = Pq 1 , recursion operator N = P\Pq and 
Darboux-Nijenhuis variables on the generic symplectic leaves of M = so* (4) for all the bivectors 
om the Section 3. 




4.2 First bivector 

Let us consider bivector Pj 
can always put 



1 ' ' " 3.3p on M — so* (4) . Using orthogonal transformations (|3.4p we 
a = (0,0, a 3 ), 6= (0,0,63). 



In this case bivector p[ X ^ gives rise to the following Darboux-Nijenhuis coordinates 

qi = 2aa 3 s 3 , q 2 = 2ab 3 t 3 , (4.6) 



and momenta 



1 



Pi = 



P2 



2a 0,3 
1 

2a 63 



arclan | — 

s 2 



f3a 3 



a u ln(a 3 s 3 - 6 3 t 3 ), 
\a l o 3 



■ arctan — 



f3a 3 



ln(a 3 s 3 - b 3 t 3 ), 



(4.7) 



t 2 J 4a 2 b 3 

which satisfy relations (|4.2p . It means that projections of Pq and Pi on the generic symplectic 
leaves of Pq are equal to 



P) 



( 


-1 

V 



10^ 
1 
000 
-1 j 



p, 



(i) 



/ 


-91 
V 



qx ^ 

q 2 



-92 / 



(4.8) 



and, therefore, symplectic foliations of Pq and Pi coincide. On the other hand it means that 
generic symplectic leaves are regular semisimple loN manifolds [3]- 
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4.3 Second bivector 

(2^ _____ _____ 

Let us consider bivector P| ' (|3.6p on M = so* (4). Using orthogonal transformations (|3.4[) and 
scaling f|3. 5|) we can always put 

a = (0,0,1), &=(6 1 ,i6 1 ,0). (4.9) 

In this case the Darboux-Nijenhuis coordinates are roots of the minimal polynomial of the 
recursion operator 

N m 

A{X) = (A — q x ){X - q 2 ) = A 2 + (t 3 - s 3 )X + (ti + if 2 )(si + zs 2 )6 2 - s 3 i 3 (4.10) 
whereas the corresponding momenta are equal to 

Pl , 2 = -i]nB(X = q 1}2 ), (4.11) 

where 

23(A) = (si - is 2 + blih + it 2 ))X + t 3 (si - is 2 ) + 6 2 s 3 (ii + it 2 ) , 

As above these variables satisfy relations (|4.2|) and, therefore, symplectic foliations of Po and 
P^ coincide. 

According [11], we can easy derive relations (14.21) from the following relations 
{A(X),B(fi)} k = ^-(A fe B(A)%)- M M(A)B( M )), 

(4.12) 

{A(X)A(n)} k = {B(X)B(n)} k = 0, fc = 0,l. 

4.4 Third bivector 

Let us consider bivector P^ (|3.7j) on M = so* (4). Using orthogonal transformations (|3.4|) and 
scaling (|3 . 5[) we can always put 

o= (0,0,1), b={b 1 ,ib 1 ,0), c= ( Cl ,-zci,0). (4.13) 

In this case the Darboux-Nijenhuis coordinates are roots of the minimal polynomial of the 
corresponding recursion operators 

A(X) = (A - qi)(X - q 2 ) = A 2 - (s 3 + i 3 )A + c x 6i(*i - it 2 )( Sl + is 2 ) + s 3 t 3 ■ (4.14) 

The corresponding momenta are defined by 

Pl , 2 = -i]nB(X = q 1}2 ), (4.15) 

where 

B(X) = ((si - is 2 ) + cibx(ti - it 2 ))X + ci&i(ii - it 2 )s 3 - (si - is 2 )i 3 . 

As above these variables satisfy relations (|4.2p and, therefore, symplectic foliations of Po and 
P x coincide. 

It is easy to see that polynomial -4(A) (|4.10p and polynomial .4(A) (|4.14[) coincide after the 
following canonical transformation 

< 2 ^-*2, t 3 ^-t 3 , at S 2 = ci6i. (4.16) 

Here b\ is the entry of the vector b (|4. 9[) , whereas b\ , cl are the entry of the vectors b and c 

(031 

So, the corresponding Darboux-Nijenhuis variables and, therefore, bivectors P x t|3. 6|) and 
P^ (|3.7p are equivalent up to canonical transformations. 
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4.5 Fourth bivector 

Let us consider bivector Q3.9[) on M = e*(3). The Darboux-Nijenhuis coordinates are roots 
of the following minimal polynomial of the recursion operator N 

A(\) = A 2 + ((a, J) - {a A b, x)aj A + - ^ - ((a, x) (b, x) - (a A xv, b A x)j 



^(a,b)(aAx,J) -^(aAx,J) + -^(x,x) -^(a,x)(b,x) (4.17) 



Linear canonical transformations of e*(3) consist of rotations 

x^aUx, J->UJ, (4.18) 
where a is an arbitrary parameter and U is an orthogonal matrix, and shifts 

x^x, J^J + Sx, (4.19) 

where S is an arbitrary 3x3 skew-symmetric constant matrix [5] . Using these canonical trans- 
formations of the vectors x and J we can always put 

A(X) = A 2 + J 3 A + T^faJi - x x J 2 ) + +x 2 2 + x\). 



However, in this case we could not directly solve equations (|4.2| with respect to the functions 

P\,2{X,J). 

Nevertheless, projections of Pq and P{ on the generic symplectic leaves of Pq are given 
by (14. 8p and, therefore, symplectic foliations of Pq and P^ coincide. 

4.6 Higher order Poisson bivectors 

For all the considered above Poisson bivectors P^ we can build the higher order Poisson 
structures by using the corresponding Darboux-Nijenhuis variables. Namely, if we restrict our 
bivectors Pq and p^™-* on their common symplectic leaves and postulate the following Poisson 
brackets between the Darboux-Nijenhuis variables 

tii, q 3 }k = {pi,Pj}k = 0, {qi,Pj} = q\hh 

then after transformations (q,p) — > (u, v) — > (s,t) one gets k+ 1 order Poisson bivectors on 
so* (4). 

As an example we present the cubic Poisson bivector 

p(3) _ f i(a,~ 1 {a,s) 2 -(b,s){c,t))s M , - (a,s+t) [a~ 1 (aAs)® (aAt) + (bAs)® (cAt)] \ (4 20) 

2 V * i(a- 1 (a,t) 2 -(f,,s)<c,i))tM / ' 

associated with the variables (|4. 14j) - (|4. 1 5[) and, therefore, compatible with the quadratic tensor 
Pi (3) (1371) . 

5 Bi-integrable systems 
5.1 The Jacobi method 

In order to get bi-integrable systems on M we can identify Darboux-Nijenhuis variables with the 
separated variables and substitute all the pairs of variables qj , pj into the separated equations 

®j(qj,Pj,ati,a 2 ) = 0, 3=1,2, (5.1) 

where $j are functions on pj, qj and two parameters a% t 2 only. 
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According to the Jacobi theorem if we solve the separated equations (15. ip with respect to 
parameters a\, 2 one gets a pair of independent integrals of motion 

Qfi,2 = Hi^(p,q), (5.2) 

as functions on the phase space M — so* (4), which arc in the bi- volution 

{H 1 ,H 2 } = {H 1 ,H 2 } 1 = 0. (5.3) 

with respect to the Poisson brackets associated with bivectors Pq and Pi (see Propositions 1-3 
in [12]). 

For instance, substituting the Darboux-Nijenhuis variables Pj,qj (|4.6H4.7p associated with 
the first bivector P^ 1 ' into the relations 

$i = 2q\ + 2V(jn) -H x - fr a , $2 = 2q\ - Hi + H 2 , 

where V = e iiaa3Pl + e - 4laa3Pl , then at /3 = -ia, b 3 = a 3 = 1 one gets 

H 1)2 = {q\ ± q 2 2 ) + V = 4a?(4 ± t\) + (i 3 - S3 )^±^ + (t 3 - S3 )" l£l -^ . 

si - «s 2 si + is 2 

Using other separated equations $1.2 = we can get many other more complicated bi-integrable 
systems on so* (4). 

5.2 Quadratic integrals of motion 

In this Section we will substitute all the known pairs of integrals of motion Hi 2 on so* (4) 
into the equations (|5.3[) and try to found bi-integrable systems associated with one of the 
Poisson bivectors from Section 3. So, in this Section we will forget about the Darboux-Nijenhuis 
coordinates and will start with integrals of motion listed in [2] and in [8] , which have different 
physical applications. 

In order to describe these bi-integrable systems we prefer to use (x, J) coordinates (|4.5[) as 
in [1(5]. 

Proposition 6 If A = —ixa, \a\ = 1 and B is an arbitrary vector, then at 

b A c + 2ia = 

the following integrals of motion 

Hi = (A, B)\J\ 2 - 2(A, J)(B, J) + (B, J Ax), 

H 2 = (B, J) (2{A, J Ax)- x 2 (J, J) + (x, x)\ 



(5.4) 



and 

Hi = {A, J A {B A J)) + (B, J A x) 



(5.5) 



H 2 = (J,B) 2 ((J AA, J A A) +2(A, J Ax) + (x,x)^ 
are in bi-involution 

{Hi,H 2 } = {Hi,H 2 }i = {Hi,H 2 } = {H 1 ,H 2 } = 0, (5.6) 
with respect to the Poisson brackets associated with Pq &2.ty) and Pi \S. 7| ). 
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Integrable system with the cubic integrals of motion H 2 (|5.4j) has been proposed in [4] . For 
this system we know the Lax matrices and the separated variables @]. The second integrable 
systems with fourth order integral of motion H 2 (|5.5[) has been considered in [7] . 

According to [3] the bi-involutivity of integrals of motion (|5.6p is equivalent to the existence 
of non-degenerate control matrices F and F, such that 



1,2. 



(5.7) 



In our case they look like 

/ (a, J) 



F = 



2k 



-iH 2 
V 2x(B,J) 







F 



/ 



2 

-iH 2 



4x(B, J) 
(a, J) 



(5.8) 



/ 



Entries of F are polynomials, whereas one of the entries of F is a rational function on so* (4). 

The eigenvalues of the matrices F and F coincide to each other. They are roots of the 
common characteristic polynomial 



A[\) = (A - Ql )(X - q 2 ) = A 2 - (a, J)A 



(J, J) i(a,x/\J) (x,x) 
~~ I 2x 4^2" 



(5.9) 



Of course, this polynomial coincides with the minimal polynomial of the recursion operator 
after suitable canonical transformation. 

Using relations (|4.12[) we can reconstruct the conjugated momenta p±_ 2 . Namely, if numer- 
ical vector d satisfies conditions (a, d) — (d, d) — 0, the coordinates q\ t2 (|5.9p and the momenta 



Pi, 2 



-<InS(A = qx <2 ), B(X) = {(d,J),A(X)} 



(5.10) 



are the Darboux-Nijenhuis variables fulfilling equations (|4.2|) . The variables p\ :2 (|5.10[) are 
defined up to canonical transformation pi — > pi + fi(gi), where fi are arbitrary functions on 
only. 

Proposition 7 Coordinates qi -2 i5.9\) and momenta p\ t2 \5.1Cfy are the separated variables for 
the integrable systems with integrals of motion H\_ 2 J5.^[ ) and H\ 2 &5.5\) . If 

b* = c, and d = c, 

the corresponding separated equations are equal to 

4k 2 (a, B) ql+q k H x -H 2 



2x 2 (c, a A B) {q\ - - C a ) q k e"** 

2x 2 {b,af\B) q k e lpk 



(5.11) 



2^(a, B) q\ 2 -HiT\ H 2 = x(c, a A B) {q 2 2 — C\)(q\ 2 — C 2 ) e~ ipi '' 

+ x(b, a A B) e ipi - 2 



(5.12) 



The separated equations (|5.11[) and (|5 . 1 2[) are related with the parabolic and cartesian Stackel 
webs on a plane. The Stackel matrices S and S diagonalize the control matrices F and F: 



F = S 



qi 
q 2 



s-\ 



s 



1 1 

2ixq 2 2ixqi 
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and 



F S(K )s-\ S=( 1 F r V 



The right hand sides of the separated equations (|5.1ip and (|5.12p are the generalized Stackel 
potentials. 

For special values of A, B and x the Hamiltonians (|5.4|) and (|5.5[) are real functions [H[7]. 
There is one integrable system with the complex Hamiltonians 



2 



H\ — a,J 2 J i + x 2 J\ — x\ J 2 , a = ix , 



a 



H 2 = a(J 1 x 2 -J 2 x 1 )(»?\J\ 2 -\x\ 2 ) + *?((J 1 x 2 - J 2 xif + {JsX! - J lX3 ) 2 ^ 
- a 2 ((J1X2 - J2X1) 2 + {J 2 x 3 - J3X2) 2 ) 

which are in the bi-involution with respect to the same Poisson brackets at a± = a 2 = and 
03 = —1. At the arbitrary value of a the Lax matrices and the separated variables for this 
system have been constructed in [§1 [3] . 

5.3 Inhomogeneous integrals of motion 

There is only one nontrivial linear Poisson bivector, which is compatible with the canonical 
bivector P an d the quadratic bivector P^ 3 ' (|3.7| 

< = (T'« 2 l). 

The linear combination 

is an inhomogeneous Poisson bivector on so* (4) compatible with Pq and such that PfdC\^ 2 = 0. 
Let us consider the following inhomogeneous Hamiltonians [4j [8] 

Hl = H x + (k,J)+a 12 (B,x), (5.14) 

Hi = H 2 + 2a 12 (B, J)(,A,x) - (k, A) J 2 - (k, J Ax) 
— K 2 a\ 2 {B, J) - ai 2 (k, x) , 

and 

Hf = Hi + {a 3 A + a A AAB,J) +a 12 {B,x), (5.15) 

HI = H 2 + 2(B 1 J)\a 3 ((A,J) 2 + x 2 \J\ 2 - \x\ 2 - 2(A, J Ax)) 

+ oti{(B A A, J Ax) + (A, J) {A, B A J)) + a 12 (B, J)(A,x) 

\ (2\x\ 2 — {A, J) 2 + 2(A, J A x)^j + 2a 4 a 12 ((B, J) (A, B A z)) 

\((B,B)(A, J) 2 - x 2 (B, J) 2 - 2(A,B)(A, J)(B, J)) 
- Aa 3 a 12 (B, J)(A,x) -a 2 12 x 2 {B, J) 2 
+ 2a 3 a 12 ((AAB, J Ax) - (A, J)(A,B A J)) 

+ 2u 3 a 12 (a 3 {A,x) + a i (A, B Ax) + a 12 x 2 {B, J) J , 

where a 12 = i(ai — a 2 ), Q!i 4 are arbitrary parameters and k is arbitrary numerical vector. 



2aj' 
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Proposition 8 The inhomogeneous integrals of motion \5. 1$ and i5.15\) are in the bi-involution 
with respect to the canonical brackets {.,.}o and the second Poisson bracket {., .}i associated 
with the bivector Pf A5.13\) 

The separated variables for these inhomogeneous bi-integrable systems are different, but they 
remain the eigenvalues of the corresponding control matrices 



F 9 = F 



and 



F g 



(a, J) + ai + a 2 



a i + a 2 
Uk, J Ax) - (k,B)\J\ 2 ai + a 2 



V 



2 ix((B,J)-a 3 ) 
iH^(a 3 — ay = 0) ({B, J)) - a 3 (a, J) + a\ + a 2 



J 



x(B, J) 2 2 

For the inhomogeneous Hamiltonians we have to add to the left hand side of the separated 
relations (|5.1ip and (|5.12p terms proportional q 2 and qk respectively. Moreover, in the right 
hand side we have to substitute ((qk — iai) 2 — Cx)((qk — ict 2 ) 2 — C 2 ) instead of (q^ — C\)(q\ — C 2 ) . 



6 Conclusion 

We classify quadratic Poisson bivectors having the given canonical foliation on M — so* (4) or 
M = e*(3) as their symplectic leaf foliation. The corresponding Darboux-Nijenhuis variables 
are constructed. A pair of known integrable systems on so* (4) can be related with one of 
such quadratic Poisson bivectors. We prove that the corresponding integrals of motion admit 
separation of variables and the separated coordinates are eigenvalues of the control matrices. 

Another approach to the construction of the quadratic and cubic Poisson bivectors P\ 
having common symplectic foliations with Pq has been proposed in [131 H4j . Among the corre- 
sponding bi-integrable systems there are generalized periodic Toda lattices, the XXX Heisen- 
berg magnet, the XXX Heisenberg magnet with boundary conditions, the Kowalevski top on 
so* (4), the Goryachev-Chaplygin gyrostat on e*(3), the Kowalevski-Chaplygin-Goryachev gy- 
rostat on e*(3) and some other well known integrable systems. Now we add to this list two 
other integrable systems on so* (4). 

The research was partially supported by the RFBR grant 06-01-00140 and grant NSc 
5403.2006.1. 
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